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Residual effects of heavy sparticles in the bottom quark Yukawa coupling: 
a comparative study for the MSSM and NMSSM 
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If the sparticles are relatively heavy (a few TeV) while the Higgs sector is not so heavy (niA is 
not so large), the Higgs boson Yukawa couplings can harbor sizable quantum effects of sparticles 
and these large residual effects may play a special role in probing supersymmetry at foreseeable 
colliders. In this work, focusing on the supersymmetric QCD effects in the hbb coupling (h is the 
lightest CP-even Higgs boson), we give a comparative study for the two popular supersymmetric 
models: the MSSM and NMSSM. While for both models the supersymmetric QCD can leave over 
large residual quantum effects in hbb coupling, the NMSSM can allow for a much broader region of 
such effects. Since these residual effects can be over 20% for the hbb coupling (and thus over 40% 
for the ratio Br(h — > bb)/Br(h — » t + t~)), future measurements may unravel the effects of heavy 
sparticles or even distinguish the two models. 

PACS numbers: 14.80.Cp,12.60.Fr,11.30.Qc 



I. INTRODUCTION 

Supersymmetry is a prime candidate for new physics 
beyond the Standard Model (SM). Among various su- 
persymmetric models the most extensively studied is the 
minimal supersymmetric model (MSSM) Another 
popular supersymmetric model, which may be equally or 
more attractive compared with the MSSM, is the next- 
to-minimal supersymmetric model (NMSSM) 0] since it 
can solve the /i-problem and alleviate the little hierar- 
chy. These models will soon be put to the test at the 
LHC. Thus, their phcnomcnological study is important 
and urgent. 

Although the most convincing evidence of supersym- 
metry is the detection of sparticle (sfermions, gauginos 
or Higgsinos) productions, the indirect probe through de- 
tecting the quantum effects of virtual sparticles in some 
measurable interactions will play a complementary role. 
If the sparticles are relatively heavy (say above a few 
TeV) and hence cannot be directly detected at the LHC, 
the indirect probe through quantum effects could be im- 
portant. For this end, the Higgs boson Yukawa interac- 
tions may play a special role since they can harbor siz- 
able quantum effects of heavy sparticles when the Higgs 
sector is not so heavy (i.e. itia is not so large). The 
dominant quantum effects of sparticles are from the su- 
persymmetric QCD interaction and in the literature the 
calculations have been performed in the MSSM for such 
supersymmetric QCD effects in the Higgs boson Yukawa 
couplings [1, IU and the associated Higgs production pro- 
cesses at the LHC The studies in the decoupling 
limit with heavy sparticles showed that for a light Higgs 
sector (mA is not large) the supersymmetric QCD can 



leave over large residual quantum effects in both the de- 
cay — s- 66 [3] and the productions at the LHC Q . Given 
the popularity of the NMSSM, it is necessary to extend 
the study to the NMSSM. This is the aim of this work. 

In this work we will focus on the supersymmetric QCD 
effects in the hbb coupling (h is the lightest CP-even neu- 
tral Higgs boson) and perform a comparative study for 
the two popular supersymmetric models: the MSSM and 
NMSSM. Such a study is interesting for two points: 

(i) The study of the NMSSM can accommodate the 
study of the MSSM, and in some limit the NMSSM 
results can reduce to the MSSM results. One can 
envisage that the supersymmetric residual effects in 
the MSSM can be magnified in the NMSSM. So the 
future measurements of such effects may be useful 
in telling the difference of the two models. 

(ii) Compared with the MSSM,_the NMSSM may pre- 
dict a different tree-level hbb coupling and different 
loop contribution. In the NMSSM the mass ma- 
trix and mixings of the Higgs bosons are enriched 
and thus the components of the lightest CP-even 
Higgs boson h are different from the MSSM val- 
ues. The residual supersymmetric QCD effects in 
the NMSSM may be larger and more interesting. 

Note that supersymmetry is a decoupling theory and all 
low-energy observables will recover their corresponding 
SM predictions when the mass scale of all supersymmet- 
ric particles (including the masses of sparticles and tua) 
take their heavy limits. The large residual quantum ef- 
fects of sparticles in the Higgs Yukawa couplings happen 
only in case that sparticles are heavy but ttia is light. If 
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both sparticles and tua take their heavy limit, the resid- 
ual effects of supersymmetry do vanish. Since so far such 
a split scenario (with light Higgs bosons and relatively 
heavy sparticles) remains possible, we should check its 
phenomenological consequence. 



II. CALCULATIONS 

We start our analysis by recapitulating the basics of 
the NMSSM. In the NMSSM a singlet Higgs superfield S 
is introduced and the Higgs terms in the superpotential 
are given by 



XSHd ■ H u — — S 3 . 



(1) 



where H u and Hd are the Higgs doublet superfields, and 
A and k are the dimensionless constants. Note that there 
is no explicit /i-term and an effective /x-parameter is gen- 
erated when the scalar component (S) of S develops a 
vev s: [i e ff — As. The corresponding soft SUSY break- 
ing terms are given by 



A x XSH d ■ H u - ^-kS 3 + h.c. 



(2) 



So the scalar Higgs potential is given by 

V F = \\H d -H u -nS 2 \ 2 + \\S\ 2 (\H d \ 2 + \H U \ 2 ) , (3) 



Vd = — (\Hd\ \H U \ — \H d ■ H u \ 



+ ^ (\H d 



\H U 



|2\ 2 



(4) 



V solt = m 2 d \H d \ 2 +m 2 u \H u \ 2 + m 2 s \S\ 2 

- (A x \SH d ■ H u + ^A K S a + ft.c.) , (5) 

where g\ and gi are the coupling constant of t/y(l) and 
SUl(2), respectively. So we can see that in the limit of 
vanishing A, k, A k and with the input of the effective 
Heff, the NMSSM can reduce to the MSSM. 

With the vevs v u , Vd and s, the scalar fields are ex- 
panded as 



H d = 

H u = 
S = 



H~d 



ifd 



■ up* 



(6) 

(7) 
(8) 



The mass eigenstates can be obtained by unitary rota- 
tions of interaction states, e.g., for the CP-even neutral 
mass eigenstates (h, Hi, H2) we have 



(h, H x , H 2 ) T = V2U H 



(9) 



So the lightest CP-even neutral Higgs boson h is com- 
posed by 



h = UgV2<t> d + UgV24> u + UgV2<j. 



(10) 



In the MSSM, h can be decomposed in the same way 
except that without the singlet component, U H is a 2 x 2 
matrix and can be parameterized in terms of a mixing 
angle a. 

The sbottom squared-mass matrix in the NMSSM is 
same as in the MSSM with fi replaced by n e ff- In the 
basis of (&l,&r), the squared-mass matrix is given by 



2 _ , M 2 m b X b 
l b \ m b X b M 2 R 



(11) 



where 



X b = A b + [i e fj tan [3 , 



Mi = 



M^+m 2 + M\{I\ - Q b s 2 w ) cos 2/3 , 



M% 



= Mi 



MlQ h s%f cos 2/3. 



(12) 



with sw = sin9w, tan/5 = v u /vd, i| and Q b being re- 
spectively the isospin and electric charge of the 6-quark, 
Mq and M b - R being the soft breaking masses, and A b 
being the soft breaking trilinear coupling. The sbottom 
mass eigenstates (61, 62) are obtained by the unitary ro- 
tation of the interaction eigenstates: 



b L 
bn 



= Zr 



— sin t 
cos 6, 



(13) 



where the unitary matrix parameterized by a mixing 
angle 0^ diagonalizes the mass matrix in Eq. (fTTj) . 

The coupling of sbottoms with the singlet component 
a of the Higgs boson h comes from the F-term of the 
superpotential, given by 



32 ™& 



-Xv u V2ab* L bji + h.c. 



(14) 



2mw cos f3 

In terms of the mass eigenstates of sbottoms, the vertex 
hbibj takes the form 



VLi 



2mw cos (5 



uM + uM + l/gc?) (15) 



where c l J is only for the NMSSM (or equivalently set 
= for the MSSM) , while and c*f are present for 
both models. They are given by 



2m b — 
2m b 



m z cos 2 /3 



(1 



m b 3 
2 ^ 2 cos 2 [3 



yli rrlj 



2m z s w 



A b Z'AZ 



li 7 2j 



16 \ ^D^D 

mi sin 2/3 



3m b 

7 2i ylj 



7 2i 7 2j 



(16) 



2m b 
m % s w 



1 - 77s 



2 - 2 sin 2/3 



7 2i 7 2j 



rvli ylj 



3m b 

+ M} (zjjz 2 J + z 2 ^ 

\f2Xmw sin/5 



<J2 



yli y2j 



7 2i 7 iy 
D ^ D , 



(17) 
(18) 
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Now we calculate the SUSY QCD corrections to the ver- 
tex hbb. At tree level it takes the same form in both the 
MSSM and NMSSM, given by 



V°- 

v hbb 



■ g2m b U(l 
2mw cos (3 



(19) 



The one-loop SUSY QCD corrections come from the 
Feynman diagrams shown in Fig[T] where we do not show 
the self-energy loop of b. In our calculations we use the 




(a) 

FIG. 1: Feynman diagrams for the SUSY QCD corrections to 
the vertex hbb. 

on-shell renormalization scheme Q and take the external 
b and b quarks on shell to get the effective vertex. With 
the corrections the effective vertex takes the form 



^ = ^T b (l + A v + A ct ), 



(20) 



where A v denotes the vertex correction from Figflja) and 
A ct is the counter term from the renormalization of mt> 
and the wave functions of b and b. They are given by 



A v - — (r l:i 4- Ul2 4- Ul3 



i y2j 

D 



x [mi (Z%Z% + Z%Z 2 ^ Cu + m- g (z x iZ 

■ r?2i r 7 li\ /-| 1 / 2 2 2 2 

+Z D Zjj \ Go I (m b ,m h ,m b ,m g ,m 



(21) 



.,_ i ■ sm(2e- b )[Bo(ml,m 2 s ,m- bi ) 



-B (m b , m g , mj 2 )] - 2m b [B[ (m b , m~ g , m- bi ) 
+B' 1 (ml,ml,ml 2 )] - 2m b m g sin(2^) 



x [B'o(ml,mg,ml 1 ) - B' (ml,m 2 g ,ml 2 )] 



(22) 



where a s is the strong coupling constant, nig is the gluino 
mass, and Bq i, B' 1 and Co,n are the scalar loop func- 
tions [8[ which can be calculated by using LoopTools 9] . 



III. NUMERICAL RESULTS 

In our calculation we use the package NMSSMTools 
[lol ] for the mass spectrum and the rotation matrix of 
the Higgs fields (to get the corresponding MSSM results 
we take the limit of very small values for A, k and A K , 
and results are checked by using the package FeynHiggs 
fill]). To study the decoupling limit of SUSY parti- 
cles, we assume all the soft breaking mass parameters 



{MQ,M^ R ,mg,A},) and the parameter /x e // are degener- 
ate, which are collectively denoted by Msusy- Then 
in the MSSM the SUSY parameters are {M S usy, M a , 
tan/3), while in the NMSSM there are three additional 
parameters (A, k, A k ). In our calculations we scan over 
these three additional parameters in the ranges 

- 0.5 < A, k < 0.5, - 500 GcV < A K < 500 GeV (23) 

Note here we did not consider a large A or k. Theo- 
retically, the requirement of perturbativity up to some 
cut-off scale will set upper bounds on A and k at weak 
scale (if the cut-off scale is chosen to be the GUT scale, 
a stringent bound A 2 + n 2 Si 0.5 is obtained Q). Phc- 
nomenologically, a large A or k will incur stringent con- 
straints from current experiments (see the last reference 
in 0). 

Before displaying the numerical results, we make some 
clarifications regarding to our numerical calculations: 

(1) In our calculations we used the package NMSSM- 
Tools [lfj which considered the loop corrections (es- 
pecially the stop/sbottom loops) to the effective po- 
tential, the masses and mixing angles of the Higgs 
bosons. These corrections are important and can- 
not be ignored. 

(2) Since we used the package NMSSMTools in which 
the Higgs mass matrices are diagonalized numer- 
ically, we did not use any approximate unitary 
transformation for our calculation. 

(3) For the b-quark mass mb, in the package NMSSM- 
Tools (and thus in our calculation), it is taken as 
the running mass rrib(Q) (we take Q = Msusy), 
which means that the sizable QCD loop effects 
are taken into account. Actually, although both 
the tree-level coupling U^U an(1 the one-loop SUSY 
QCD contributions SV hb b are proportional to nib 
and thus very sensitive to the value of nib, the rela- 
tive correction effects 3V hb i/V? bb - (displayed in our 
numerical results) are not so sensitive to the value 
of nib- 

Now we present some numerical results. To see the gen- 
eral feature of the corrections, we first switch off the ex- 
perimental constraints on the parameter space and in 
Figs (21 and display the resluts of the relative correction 
effects 



A 



SQCD 



~ V°- 

hbb 



(24) 



We see that for a light niA the SUSY QCD with a large 
Msusy can leave over sizable effects in the coupling hbb. 
But as niA gets heavy, such residual effects of SUSY QCD 
become small in magnitude, showing the decoupling be- 
havior of supersymmetry. Compared with the MSSM 
results, the effects in the NMSSM can vary in a much 
broader region. While the corrections are always nega- 
tive in the MSSM, in the NMSSM the corrections can be 
both negative and positive. 
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FIG. 2: SUSY QCD corrections to hbb coupling versus rriA in 
the MSSM and NMSSM for tan/3 = 20. The scatter plots are 
for the NMSSM while the curves for the MSSM. 
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FIG. 3: Same as Fig. O but versus Msusy for fixed values 
of tua- 



In order to figure out in what areas of the NMSSM pa- 
rameter space the one-loop SUSY QCD effects on the hbb 
coupling are sizable, we present a set of sample points in 
Table Q] We see that for the SUSY QCD loop effects on 
the hbb coupling to be sizable in the NMSSM, the ratio 
U^/UR ( UR 2 and llR are respectively the components 
of 4>u an d (\>d in h, as defined in EcillOp plays the key role 
for the following reasons. As shown in Eq.(|24p. Asqcd 
is composed of two parts: the counter-term part A ct and 
the vertex-loop part A v . We found that in most of the 
parameter space allowed by the LEP constraints, A ct is 
negative and dominant in size, which is independent of 
how h is composed of. Whereas, A v is positive and can- 
cel A ct to some extent. The size of A v can be enhanced 
by the ratio U^/UR, as shown in Eg . (f2"Tj) . Although the 



TABLE I: Some sample points in the parameter space of the 
NMSSM, in which the SUSY QCD loop effects Asqcd can 
be of very different size. The mass parameters are in unit of 
GeV. Here we fixed M A = 200 GeV and tan/3 = 20. With 
such fixed parameters and Msusy varying from 400 GeV to 
2 TeV, the value of m h in the MSSM varies from 106 GeV 
to 119 GeV while the corresponding values of Asqcd in the 
MSSM are about -11%. 



A 


K 


A K 


Msusy 


uli/uR 


UR/UR 


m h 


Asqcd 
(%) 


0.106 


-0.215 


446 


1788 


16.2 


-0.023 


117.4 


-4.1 


-0.116 


0.180 


333 


1575 


14.7 


0.042 


116.3 


-6.0 


0.133 


-0.256 


393 


1145 


14.8 


-0.046 


114.7 


-6.1 


0.123 


-0.164 


460 


1568 


14.1 


-0.062 


116.1 


-6.8 


-0.103 


-0.072 


-294 


1544 


10.1 


0.118 


114.8 


-12.0 


0.107 


0.073 


-175 


1831 


9.7 


-0.100 


115.2 


-12.4 


0.126 


0.101 


-426 


1738 


7.3 


-0.069 


114.8 


-15.5 


-0.106 


-0.123 


-377 


1401 


5.4 


0.024 


115.3 


-18.0 


0.110 


0.130 


-324 


1419 


4.1 


-0.017 


114.9 


-19.7 


0.106 


0.406 


-304 


623 


0.5 


0.002 


99.4 


-24.5 



ratio U^/UR can also enhance the size of A v , its effect 
is suppressed by the smallness of and A in c l J (in our 
scan we found that both A and are small in order 
to satisfy the LEP constraints encoded in the NMSSM- 
Tools). Therefore, as U[i,/Uf{ gets large, A v becomes 
more sizable and, due to its cancellation effect, the total 
correction effects become less sizable. Note that although 
the value oirrih is below 114 GeV for some points, these 
points can still satisfy the LEP constraints encoded in 
the NMSSMTools because the LEP bound on the MSSM 
mh is 92 GeV (see Fig. 2 in the paper by Barger et. al. in 
Ref . Q ) • By the way, some comprehensive studies on the 
phenomenology of the NMSSM Higgs sector (which may 
be quite different from the phenomenology of the MSSM 
Higgs sector) have been performed in Ref.JH, where all 
the Higgs couplings and all the decay modes as well as 
productions have been intensively studied. 

From Table U we see that although A is scanned in the 
range of —0.5 ~ 0.5, only small values of A (~ 0.1 in 
magnitude) survived the LEP constraints encoded in the 
NMSSMTools. As A gets large, the mixing between sin- 
glet and doublet Higgs fields becomes large (although the 
mixing is not totally determined by A and other param- 
eters are also contributing) and, consequently, is more 
constrained by LEP experiments 0). Such a small A, to- 
gether with other NMSSM parameters appearing in the 
mass matrix of the Higgs fields, leads to a small mixing 
between singlet and doublet Higgs fields. As shown in 
Table HI the singlet component in h is small, which 
means that the lightest Higgs boson h is not the singlet- 
like one and instead it is doublet-dominant (MSSM-like). 
For the hbb coupling with such a MSSM-like h, the 
NMSSM can still allow for quite different SUSY QCD 
effects compared with the MSSM predictions. The rea- 
son is that, as discussed above, the SUSY QCD effects 
are sensitive to the ratio U^/UH, which is different in 
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these two models (even for a small A, the NMSSM can 
still allow for different-size components and Uf{ in h 
because other NMSSM parameters are also contributing 
to the mass matrix and mixings of the Higgs fields). 

When we switch on the comprehensive experimental 
constraints encoded in the NMSSMTools [l(|, we obtain 
the results of Asqcd displayed in FigJU These experi- 
mental constraints are comprehensive, including the LEP 
II searches for the Higgs boson, various B-decays and the 
muon anomalous magnetic moment (muon g — 2). 




1800 2000 



SUSY 



(GeV) 



FIG. 4: Same as Fig. El but for m A = 100 GeV in the 
NMSSM. The dark-shaded region denoted by 'x' (red in 
color) are allowed by the experimental constraints considered 
in the NMSSMTools [hJ. The light-shaded region denoted 
by '•' (yellow in color) are excluded by LEP experiment or B 
physics, and the medium-shaded region denoted by '+' (sky- 
blue in color) are excluded by the muon g — 2. 

From Fig|4] we see that in the special scenario under 
our consideration (all soft breaking mass parameters are 
degenerate), Msusy is constrained in a certain range. 
The range of Msusy 600 GeV is not allowed by LEP 
experiments, while Msusy 1 TeV cannot explain the 
muon g — 2 data (we require the supersymmetric effects 
to account for the deviation of a e * p — a^ M = (29.5 ± 
8.8) x 10~ 10 at 2er level). In the allowed region of the 
parameter space, the SUSY QCD corrections to the hbb 
coupling can still be significant, over 20% in magnitude. 
If we switch off the muon g — 2 constraint ( the hadronic 



contribution to 



,SAf 



is not so certain [12|]), the allowed 



parameter space gets much broader, as shown in Figf4] In 
our study we did not require supersymmetry to explain 
various plausible evidences of dark matter. 

Note that the SUSY contributions to the muon g— 2 are 
sensitive to the soft masses in slepton sector and chargino 
sector, but not dependent on the squark or gluino mass 
involved in our SUSY QCD loops. The stringent con- 



straint from the muon g — 2 data shown in Figfi] comes 
from our simple assumption that all soft masses (in the 
squark sector, the slepton sector and the gaugino sector) 
are degenerate. Of course, the muon g — 2 constraint is 
not necessary; without such an assumption of degeneracy, 
the constraint is lifted. 

Since the SUSY QCD residual effects can be over 20% 
for the hbb coupling and thus over 40% for the ratio of 
the branching fractions (SUSY QCD does not contribute 
to h — > t + t~ at one-loop level) 



Rb/r — 



Br(h -> bb) 
Br(h -> t+t- 



(25) 



future measurements at the LHC or ILC may unravel 
such supersymmetric effects. This ratio Rb/ T is proposed 
in [f| as a probe to new physics. To measure this ratio 
at the LHC, one may count the event numbers of the pro- 
duction of hbb followed respectively by the decay h — ► bb 
and h — ► t + t~, and the ratio of these two event-numbers 
can be a measure of Rb/ T (the difference of efficiency for 
6-tagging and r-tagging should be taken into account). 

Finally, we make some remarks regarding to our re- 
sults: 

(1) We only investigated the SUSY QCD corrections, 
which is 0(a s ) and thus should be the most im- 
portant among the SUSY corrections. Among 
the SUSY electroweak corrections, the Higgs-top 

2 

Yukawa corrections, which are Ota— may also 

be sizable although they are seemingly not as large 
as the 0(a s ) SUSY QCD corrections. With the 
consideration of the SUSY electroweak corrections, 
our conclusion will remain qualatively unchanged. 

(2) Since our main interest is the residual effects 
of heavy sparticles, we assumed all soft SUSY- 
breaking mass parameters are equal to Msusy- 
This is a very strong simplifying assumption. If 
we lift such an assumption and consider the mul- 
tiple free soft parameters, we will obtain the re- 
sults which numerically are different to some ex- 
tent while qualitatively exhibit the same feature, 
i.e., the NMSSM allows broader corrections than 
the MSSM, because the NMSSM parameters can 
complicate the mass matrix and mixings of Higgs 
fields. 

(3) We should stress again that supersymmetry is a 
decoupling theory and the large residual quantum 
effects of sparticles in the Higgs Yukawa couplings 
are present only in case of a light m^. If both 
Msusy and itia take their heavy limit, the SUSY 
effects will vanish. Such a decoupling behavior is 
similar in the NMSSM and MSSM, as shown in 
Fig.2. 
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IV. SUMMARY 

We focused on the SUSY QCD effects in the hbb cou- 
pling and performed a comparative study for the two 
popular SUSY models: the MSSM and NMSSM. We 
found that for both models the SUSY QCD can leave 
over large residual quantum effects in hbb coupling if 
the sparticles are relatively heavy (a few TeV) while the 
Higgs sector is not so heavy (m^ is not so large). Com- 
pared with the MSSM results, the NMSSM can allow 
for a much broader region of such residual effects. Since 
these residual effects can be over 20% in magnitude, fu- 



ture measurements of the hbb coupling may unravel such 
supersymmctric effects or even distinguish the two mod- 
els. 
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